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ABSTRACT 

Gamma Ray Bursts (GRBs) show evidence of different spectral shapes, hght curves, duration, 
host galaxies and they explode within a wide redshift range. However, the most of them seems 
to follow very tight correlations among some observed quantities relating to their energetic. 
If true, these correlations have significant implications on burst physics, giving constraints on 
theoretical models. Moreover, several suggestions have been made to use these correlations in 
order to calibrate GRBs as standard candles and to constrain the cosmological parameters. 
We investigate the cosmological relation between low energy a index in GRBs prompt spectra 
and the redshift z- We present a statistical analysis of the relation between the total isotropic 
energy Ei,,, and the peak energy Ej, (also known as Amati relation) in GRBs spectra searching 
for possible functional biases. 

Possible implications on the Etso vs Ep relation of the a vs (1 +z) correlation are evaluated. We 
used MonteCarlo simulations and the boostrap method to evaluate how large are the effects of 
functional biases on the E/so vs Ep. We show that high values of the linear correlation coeffi- 
cent, up to about 0.8, in the E/so vs Ep relation are obtained for random generated samples of 
GRBs, confirming the relevance of functional biases. 

Astrophysical consequences from Ej^o vs Ep relation are then to be revised after a more accu- 
rate and possibly bias free analysis. 

Key words: gamma-rays: bursts, gamma-rays: observations. X-rays: general, methods: sta- 
tistical. 



1 INTRODUCTION 

Gamma-ray Bursts (GRBs) are brief and intense flashes of high en- 
ergy radiation emitted mostly in the y-ray band. They are detected 
from wholly random directions in the sky at the rate of about once 
a day and typically last from a few milliseconds to several minutes. 
Within a few years, the BATSE experiment on board the NASA's 
Compton Gamma Ray Observatory satellite (Fishman et al. 1989) 
has recorded over 2700 GRB events with an isotropic distribution 
in the sky (Meegan et al. 1996). However, although BATSE was 
very sensitive to high-energy photons, it could not discern the di- 
rection of a burst to better than a few degrees uncertainty, too large 
to pinpoint the location of individual explosions. 

The real revolutionary step forward occurred in the 1997, 
thanks to the Itahan-Dutch BeppoSAX satellite (Boella et al. 1997). 
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This satellite was not as sensitive as BATSE to y rays, but its rela- 
tive quick response of pointing system, coupled with good accuracy 
position information, permitted the first detection of an X-ray af- 
terglow, the radiation emitted after the initial burst of y-ray (Costa 
et al. 1997). This discovery of afterglows made redshift measure- 
ments possible and confirmed that GRBs lie at cosmological dis- 
tance (0.0085 (Galama et al. 1999) <z< 6.29 (Kawai et al. 2006)). 

It is well known that in the past years several correlations have 
been discovered linking various energies characterizing GRBs. All 
of them involve Ep, the energy peak of time integrated spectral en- 
ergy distribution (SED). The first relation found links the rest frame 
isotropic energy £,„ with E,, (Amati et al. 2002, AM02 hereinafter, 
also known as 'Amati relation', see also Amati, 2006, AM06 here- 
inafter, for an updated version). This correlation is seen by several 
authors as an useful method to standardizing GRB energetics. 

Subsequently, Ep was found to be tightly correlated also with 
the collimation corrected energy Ey (Ghirlanda et al. 2004, GH04 
hereinafter, also known as 'Ghirlanda relation'), and this relation is 
used to constrain cosmological parameters using GRBs as 'known' 
candles. The same relations can be expressed in terms of luminosity 
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using spectra not time integrated: L,,o oc (Yonetoku et al. 2004), 
LyOcE'l, (Ghirlanda, Ghisellini & Firmani 2006). 

However, this kind of relations seems to contradict some ob- 
servational evidences, as the large variety of light curves, spectra, 
redshifts, durations and host galaxies, leading us to suppose that the 
nature of GRB explosions is not unique. 

Band et al. (2005) and Nakar et al. (2005) tested the consis- 
tency of a large sample of BATSE GRBs (with unknown redshift) 
with the Amati and Ghirlanda relations. Their results suggest that 
they may be artifacts of selection effects, inferring that about half 
(Nakar et al. 2005) or even ~ 90% (Band et al. 2005) of the whole 
GRB population cannot satisfy the correlation for any value of the 
redshift. However, these conclusions have been questioned by sev- 
eral other authors (Ghirlanda et al. 2005, Bosnjak et al. 2005, Pizzi- 
chini et al. 2005), that found instead that the peak energy and the 
fluences of BATSE GRBs with unknown redshift are fully consis- 
tent with the Eiso vs Ep and the Ey vs Ep relations. 

In AM02, and in AM06, Eiio and Ep are not evaluated directly 
by a fitting procedure but they are calculated using analytic rela- 
tions that include the same parameters. These procedures could in- 
troduce functional biases in the spectral correlation quoted above. 
The influence of these functional biases has not been never consid- 
ered in literature. 

In this work we present a statistical analysis of the £,,„ vs Ep 
relation, based on Monte Carlo and bootstrap simulations, search- 
ing for possible intrinsic correlation terms. We previously study 
the cosmological relation between the redshift z and the low energy 
spectral index a (see Sec. 3), because a correlation between these 
spectral parameters can make more tight the vs Ep relation. In 
Sect. 4 we evaluate functional and intrinsic correlation terms in the 
Eiso vs Ep relation, in order to estimate if the method used to derive 
it is statistically biased. 



energy Ep of the SED is related to the spectral parameters by: 

Ep = (a + 2)Ec < Ec (4) 
being, typically, -2 < a < -1. 



3 THE CORRELATION BETWEEN a AND THE 
REDSHIFT 

For the GRBs of BATSE catalogue a correlation between a and 
the redshift z was found (Lloyd et al. 2000). It was interpreted 
as partially due to the dependence of the a on the spectral cur- 
vature around Ep (Lloyd et al. 2000) and partially to the depen- 
dence of Ep on the a index (Band et al. 1993). A similar correlation 
was found in AM02 with a high logarithmic correlation coefficient 
''log = -0.83 between the photon index a and (1 + z) (which cor- 
responds to a linear correlation coefficient of the same quantities 
riin = -0.77). No clear explanation was given by the authors, who 
suggested that it could be a consequence of the Ei^g vs Ep corre- 
lation. Even if a real satisfactory explanation for the a vs (1 -I- z) 
relation has not yet been found, we will show (Sect. 4), that it is 
relevant when studying the Amati and similar relations. 

We first searched for its validity using the GRB samples in 
GH04 and in AM06, that is larger than the one of AM02. After 
excluding the bursts of AM02, because their parameters were eval- 
uated in the GRB rest frame, we obtained two samples of 15 and 23 
bursts, respectively. For both samples we found a significant corre- 
lation, with rij„ = -0.61 for GH04 bursts and r;,„ = -0.56 for the 
AM06 sample, having the same trend of AM02. Logharitm correla- 
tion coefficents and p-chance probabihty are reported in Tab.l. We 
fitted the data in AM02, GH04 and AM06 samples with a linear 
relation: 



■ m(l + z) + q 



(5) 



2 THE GRB SPECTRAL DESCRIPTION 

GRBs have a non thermal spectrum that varies strongly from one 
burst to another. It is generally found that a simple power law does 
not fit well their spectra because of a continuos steepening toward 
the high energies. An excellent phenomenological model was intro- 
duced by Band et al. (1993) to describe the prompt time-integrated 
GRBs spectra, composed by two power laws joined smoothly at a 
break energy Ei,: 



N{E): 



E<E, 
E>Eh 



(1) 



where N{E) is the number of photons per unit of area and energy, 
while Eq is a reference energy usually fixed to the value of 100 keV. 
Under the continuity requirement for the function N{E) and its first 
derivative, the break energy and normaUzation are given by: 



E,, = {a-P)Ec 
for typical values a ^ -1.4 and p ^ -2.4, Ei, » Ec, and 

-la-p 



B = A 



ia-P)Ec 



(2) 



(3) 



There is no particular theoretical model that predicts this spec- 
tral shape, however it provides good fits to most of the observed 
spectra in terms of four parameters: the two photon indices a and 
j8, the exponential cut-off E^ and the normaUzation constant A, pa- 
rameters directly estimated during the fitting procedure. The peak 



The best fit values of parameters m and q are reported in Tab.l, 
while regression lines are plotted in Fig.l for each GRB samples. 

A possible explanation of the a vs (1 -i- z) correlation could be 
due to a selection effect: brighest GRBs correspond to flatter photon 
indicies a. The correlation between a and the redshift implies that 
we are loosing faint GRBs at high z values. In Fig. 1 the GRBs 
number decreasing at high redshift is shown. 
In the next section we show how this correlation affects the £,.o vs 



4 THE E,so VS Ep RELATION 

It is clearly known that when two or more uncorrelated and inde- 
pendent parameters are used to build up other variables, a correla- 
tion term will arise between these new variables due to the func- 
tional relations used in the calculations. So, in searching any kind 
of correlations between these new variables it is necessary to take 
into account of the possible biases due to the method adopted to 
estimated them. 

AM02 reported the analysis of 12 GRBs with known redshifts 
observed by BeppoSAX. For 9 blue-shifted spectra, using the Band 
model, the spectral parameters {a, fl, E^ and A) are estimated in 
each GRB rest frame. The authors found a very tight correlation 
(r;i„ = 0.96) between Ep, defined by Eq.4 and: 



(l+z)2 



/ 



ifkeV 



ENiE;a,/3,Ec,A)dE 



(6) 



Functional biases in GRB 's spectral parameter correlations 3 



Table 1. a-x correlation 



Sample REF. GRBs(#) z range r;,„ P;,„ riog Pi„g q m 



Amati et al., 2002 9 0.42- 3.42 0.77 (0.10) 0.006 -0.83 (0.19) 0.002 -2.0688 0.32312 

Ghirlandaetal.,2004 15 0.0085-4.5 0.61(0.05) 0.01 -0.67(0.15) 0.004 -1.466 0.14982 

Amati, 2006 23 0.0085-4.5 0.56(0.04) 0.04 -0.63 (0.05) 0.0008 -1.3813 0.21449 



I AMATl 2002 | 




GHIRLANDA 2004 
AMATI 2006 



(1+z) 



Figure 1. The vs (1 + z) relation. In the upper panel the blueshifted GRBs 
(black filled circles) of the AM02 sample are reported, while, in the lower 
panel, there are GH04(black filled circles) and AM06(red crosses) bursts. In 
both panels the regression lines are plotted. 



the total (isotropic) energy emitted in the GRB rest frame. Both 
quantities depend on the two parameter a and Ec and it is worth to 
be noted that is a function of z, that is related to a as discussed 
in the previous section. 

In our analysis a flat Friedmann-Robertson- Walker Universe 
is assumed, with Hq = 65 km/(s Mpc), = 0.3 and = 0.7 as 
reported in AM02. The luminosity distance was calculate using 
the formula (Hogg 1999, Carroll et al. 1992): 



D, 



(1+z) 



- r 

^^0 Jo 



dz 



(7) 



that was also successfully compared with other cosmological calcu- 
lators (N. Wright and R. Priddey). Our code was previously tested 
using as input data the parameters reported in Tab. 2 of AM02 and 
finding the same results for £,s„ and 

In a recent work (Amati et al. 2006 and references therein), 
the vs correlation is evaluated using the following relations 
that do not require bluesfhited GRBs spectra: 



(8) 



We also tested functional biases in the method performed with Eq.8 
and Eq.9. 



4.1 Monte Carlo simulations: the £■,„, vs Ep correlation 

Our investigation of a possible functional bias was made by build- 
ing up a numerical code that, after assigned the five spectral param- 
eter z, a, yS, Ec, A, calculates and Ep using Eq.4 and Eq.7. 

To clarify the efi'ect of redshift on the correlation, we initially 
fixed its value at z = 1, while the other spectral parameters a, f3, E^, 
A were randomly, and so uncorrelated, generated with an uniform 
distributions within the intervals of AM02 [-0.7; -1.9], [-2.1; -2.7], 
[340; 840] keV, [0.07; 2.27] ph/Ccm^keV), respectively. Different 
values of z were tested, observing that the correlation do not de- 
pend on its choice. Our analysis was also performed with ranges of 
spectral parameters wider than the observed ones, with consistent 
results. Ej,o and Ep values were calculated for each set of spectral 
parameters, iterating this procedure for 45 GRBs, a factor of five 
larger than those analysed in AM02. Values of Ep < 90 keV were 
discarded to be in agreement with AM02. 

Finally, the correlation coefficients r;,„, ri„g and the non- 
parametric Spearman correlation coefficent f jp, , with their standard 
deviations were estimated for this sample. This procedure was re- 
peated 1000 times to study the distribution of the correlation coeffi- 
cients. The entire procedure have been also tested using 100 GRBs 
in each samle with consistent results with the previous analysis. 

In Fig.2 (upper panel) the distribution of r;,„ is represented, 
with a mean value of < r;;„ >= 0.48 and a standard deviation of 
o";,„ = 0.13. The results obtained using the logarithmic coefficient 
have a mean value < ri^g >= 0.26 and cr/og = 0.14, while the mean 
Spearman correlation coefficent is < r,pr >= 0.33 with a standard 
deviation cr^pr = 0.10. The mean chance probability for the linear 
correlation coefficent is 0.03 and the probability to find a spuri- 
ous correlation coefficent is not negligible, being a relatively large 
number of linear correlation coefficents higher than 0.6. 

Sakamoto et al. (2006) shows for a sample of 32 GRBs that 
there is a correlation between the energy peak and the total fluence 
with riog = 0.58 Our simulations at fixed redshift, corresponding 
to the relation between these quantities, show that we found a cor- 
relation coefficient practically coincident with the Sakamoto et al. 
results. Again this support the relevance of biases introduced by 
functional relation. 

To evaluate the cosmological influence due to the redshift on 
the Eiso vs Ep correlation, we also generated the z values randomly. 
To taking into account the correlation relating to a and z, as dis- 
cussed in the previous section, redshift values were generated with 
an uniform distribution for each GRB in the range [0.42; 3.42], 
while a indices were evaluated using the formula: 



I 

J Ik 



lO*l<eV/(l+z) 



(1 + Z) JikeV/{l+z) 



EN(E;a,p,E„A)dE. 



(9) 



a = m(l+z) + q + R (10) 
The values of m and q are evaluated using AM02 data (see values 
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Figure 2. The distribution of 1000 linear correlation coefficents evaluated 
with fixed redshift (upper panel) or taking into account of the intrinsic cor- 
relation between a and z (lower panel) using Eq.4 and Eq.5 for the Monte 
carlo simulations. 



Figure 3. The distribution of 1000 linear correlation coefficents evaluated 
using Eq.8 and Eq.9 instead of Eq.4 and Eq.5 and taking into account of 
the intrinsic correlation between a and z in GH04 and AM06 for the Monte 
carlo simulations. 



reported in Tab. 1), while is a random number within the range 
[-0.5,0.5]. The introduction of R allows us to have a correlation 
coefficient < r(i„ >= -0.70 between a and (1 +z) consistent with the 
observed value (AM02). The p, E^, A parameters were generated 
as in the previous analysis. 

By iterating this procedure to build 1000 correlation coeffi- 
cients, we found a value of < r;i„ >= 0.66 with o"/,„ = 0.09 
(Fig.2, lower panel) with a mean value of the chance probability 
of 3.7 X lO"'* and a value of < nog >= 0.54 with a-iog = 0.12 and 
< Tjpr >= 0.55 and cr,,,,- = 0.10. 

Following GH04 and AM06 relations taking into account of 
the intrinsic correlation between a and (1 + z) and of the different 
range of redshift (see Tab. 1), we have found an even more higher 
correlation coefficents, as reported in Tab.2. This is beacuse Ei^o 
and Ep in Eq.8 and Eq.9 depend both on the a index, the high en- 
ergy cut-off £<; (Eq.5), and by the redshift z. The distributions of 
linear correlation coefficents are reported in Fig.3. 

It is worth to be noted that if you decrease the number of sim- 
ulated GRBs to build a single correlation coefficient, for example 
from 45 down to 9 as the AM02 sample, the distribution of be- 
comes more spread and the frequency of values close to 1 is not 
negligible (r/,„ > 0.7 for about half values). 

In Tab.2 all the correlation coefficents values, calculated tak- 
ing into account of the a vs (1 -l- z) correlation, and their p-chance 
probabilities are summarized. 

In the correlations analysis could be relevant the ranges of the 
quantities involved in the evaluation of spectral parameters. To ver- 
ify our results we performed the same simulations using parameter 
ranges wider than the observed ones, for example in the case of 
the normalization we used an interval spanning up to five order of 
magnitude, and we obtain the same results. 



Table 2. Correlation coefficents 



Relations 


run 


0"K/i 


Plin 


nog 


O-log 


^spr 


O'spr 


AM02 


0.66 


0.09 


3.7 X lO-'* 


0.54 


0.12 


0.55 


0.10 


GH04 


0.58 


0.11 


4.7 X 10-3 


0.51 


0.12 


0.37 


0.13 


AM06 


0.55 


0.12 


7.9 X 10-3 


0.48 


0.12 


0.35 


0.13 



4.2 Monte Carlo simulations: the k exponent in the Ep oc 
relation 

We also checked our results evaluating the exponnt in the Eis„ vs 
Ep relation. The k exponent in the Ep oc E'l^^ relation evaluated in 
AM02 with 9 GRBs has a value of 0.52. 

In Fig. 4 we report our distribution of the k exponents in the 
relation cc with a mean value of < A: >= 0.26 and a disper- 
sion of CTj. = 0.05. 

We observe that in AM02 GRBs sample there are two pecu- 
liar GRBs: GRB 990123 having Ep outside of the observed energy 
range and GRB 010222 for which there is only a lower limit on 
Ep based on an assumption on the p index. Excluding from the 
AM02 sample these two GRBs, the obtained value for the k expo- 
nent is 0.36, in agreement with the analysis performed by Amati et 
al. (2003, 2006), and it is consistent with our results. It is worth to 
be noted also that values of k exponents in Amati et al. (2006) are 
in the range [0.35; 0.57] and seem to depend by the sample used to 
evaluate the Ej^o vs Ep relation. 
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Figure 4. The distribution of k exponents in tlie simulated GRBs of the Ejjo 
vs Ep relation. 

4.3 Boostrap method 

Our results were also checked performing also a bootstrap method, 
a techinique frequently adopted in the statistical analysis of corre- 
lations (Manly, 1997), that consists in the random association of the 
parameter's values between the elements in the same sample. 

To perform this method we randomly selected the spectral pa- 
rameters z, yS, Ec, A of eight GRBs of the sample used in AM02 
(see Tab.l and Tab.2 in AM02). Wc excluded GRBO 10222 because 
in AM02 only a lower limit of the value of the Ep is reported. We 
found values of the logarithm correlation coefflcent < >= 0.23 
with (Tiog = 0.44. 

To take into account the correlation between a and z, their val- 
ues were selected belong to the same ORB. Following this method, 
we found an higher value of < riog >= 0.52 with a standard devia- 
tion of CTiog = 0.24. 

The results of the bootstrap method (Fig. 5) are consistent with 
the Monte Carlo simulations. In Fig. 6 the coverage of our bootstrap 
simulated GRBs in comparison with observed ones (AM02) in the 
Eiso vs Ep plane is shown. 



5 DISCUSSION 

In this paper a statistical analysis of the GRBs parameters corre- 
lations a vs (1 -I- z) and Ej^o vs Ep is presented. We studied the a 
vs (1 + z) correlation observing that it is present not only in the 
AM02 GRB sample but also in those analysed by GH04. As stated 
by different authors (Lloyd et al. 2000, Band et al. 1993), there is 
no clear explanation for it, but we suggested that it could be a direct 
observational consequence of a selection effects: faint high redshift 
GRBs are under the detectability threshold. As a consequence of 
this, the decrease of high redshift GRBs introduce a correlation be- 
tween a and z. This correlation can have important consequences 
when we attempt to use GRBs as standard candles and, moreover, it 



Figure 5. The distribution of 1000 logarithm correlation coefficents for the 
boostrap simulations (upper panel) and taking into account of the intrinsic 
correlation between a and z (lower panel). 



introduces a bias in searching relations between quantities builded 

using a and z- 

We investigated how the presence of functional correlation 
terms affects the Ei^o vs Ep relation. As builded in AM02, £,jo and 
Ep depend on a common set of spectral parameters. We performed 
several Monte Carlo simulations , initially considering only the in- 
trinsic correlation term at a fixed redshift. We found that the 
vs Ep relation is biased, with a not negligible probability to repro- 
duce correlation coefficents higher than 0.5. When we introduced 
in the generation of simulated sample also the or vs (1 -I- z) corre- 
lation, the mean value of the correlation coefficent increase up to 
0.64. Higher values were obtained using Eq.8 and Eq.9 instead of 
Eq.4 and Eq.5. A boostrap method was also applied to the AM02 
sample. The good agreement with the previous analysis proves that 
our results are independent on the simulation code. 

An already discussed (e.g. Band et al. 2005) relevant subject 
is the contribution from selection effects in making more tight the 
Eiso vs Ep relation with respect to our simulations. It is possible 
that instrumental sensivity limits would produce forbidden regions 
in the plane E,™, E,,. As an example, a combination of the spectral 
parameters in the GRB rest frame could make a GRB with a very 
low luminosity and therefore below the detectability threshold. 
We also stress that the SED energy peak is an averaged quantity, 
because during a single GRB that shows more bumps in its 
lightcurve we expect that Ep changes, making difficult the link of 
this variable to physical quantities. 

In the literature appear several other tight correlations con- 
cerning the energetics and the spectral distributions of GRBs (e.g 
E^ vs Ep (GH04), L vs Ep (Yonetoku et al. 2004)), some of them 
are used to constrain cosmological parameters, using GRBs as 
'known' candles. All of them are based on the Ei^o vs Ep relation 
and this implies that they are affected by functional biases. In the 
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Figure 6. Black filled points represent 8 GRBs simulated witli the boost- 
rap method, in comparison with those in AM02 sample indicated with red 
circles. 



GH04 relation, for example, the collimation corrected energy is 
computed considering the jet opening angle, that is in turn derived 
from Ejso and z- 

In any case, we do not exclude that a physical relation between 
the Eiso and E,, can really exist, but it can be safely established only 
after removing all the biases. 

A good solution could be to use spectral laws explicitly written 
in terms of the quantities for which a correlation is searched (as 
made, for example, in Tramecere, Massaro & Cavaliere (2007) for 
the spectrum of the near HBL object Mkn 421). 

A further step could also be to work with homogeneous sam- 
ples of GRBs in terms of some relevant characteristics (e.g. time 
evolution, spectral shape, ...) to verify that the correlation is actu- 
ally followed by them. 
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